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Abstract. Dynamics of a classical particle in a one-dimensional, randomly driven 
potential is analysed both analytically and numerically. The potential considered 
here is composed of two identical spatially-periodic saw-tooth-like components, one of 
which is externally driven by a random force. We show that under certain conditions 
the particle may travel against the averaged external force performing a saltatory 
unidirectional drift with a constant velocity. Such a behavior persists also in situations 
when the external force averages out to zero. We demonstrate that the physics behind 
this phenomenon stems from a particular behavior of fluctuations in random force: 
upon reaching a certain level, random fluctuations exercise a locking function creating 
points of irreversibility which the particle can not overpass. Repeated (randomly) in 
each cycle, this results in a saltatory unidirectional drift. This mechanism resembles 
the work of an escapement-type device in watches. Considering the overdamped limit, 
we propose simple analytical estimates for the particle's terminal velocity. 
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1. Introduction 

A question of liow to make use of unbiased random or periodic fluctuations rectifying 
them into directional motion has been a long-standing challenge for mechanists, 
physicists and biologists. For fluctuations acting on the microscopic level, originating, 
e.g., from Brownian noise, this problem has been under debate since the early works 
of Maxwell and Smoluchowski |2] , who formulated several important concepts (see 
also Ref. 0). In particular, it has been realized already by Curie more than hundred 
years ago ^:\, that although the violation of the x —>■ —x symmetry is not sufficient to 
cause a net directional transport of a particle subject to a spatially asymmetric, but 
on large scale homogeneous, potential, the additional breaking of time reversal t — > — t 
symmetry, (e.g., due to dissipation), may give rise to a macroscopic net velocity. Thus, 
in principle, on the microscopic scale directed drift motion may emerge in the absence 
of any external net force. These systems, referred nowadays to as thermal ratchets [2], 
have been scrutinized both theoretically (see, e.g. Refs. [5-16]) and experimentally (see, 
e.g. Refs. [17-20]). Considerable progress in this field, useful concepts and important 
results have been summarized in several extensive reviews [21-24]. 

As a matter of fact, on the macroscopic scale of our everyday life this very problem 
of obtaining useful work from random or periodic perturbations has been and continues 
to be successfully accomplished by engineering some clever technical devices; water- 
or windmills and watches being just a few stray examples. In particular, to make 
the watches work, the watchmakers had to invent a device capable to convert the raw 
power of the driving force into regular and uniform impulses, which was realized by 
creating various so-called escapement devices (see, e.g., an exposition in the web-site 
|25j). In practice, this is most often a sort of a shaft or an arm carrying two tongues 
- pallets, which alternately engage with the teeth of a crown-wheel. The pallets follow 
the oscillating motion of the controller - the balance-wheel or a pendulum, and in each 
cycle of the controller the crown-wheel turns freely only when both pallets are out of 
contact with it. Upon contacts, the pallets provide impulses to the crown-wheel, (which 
is necessary to keep the controller from drifting to a halt), and moreover, perform a 
locking function stopping the train of wheels until the swing of the controller brings 
round the next period of release. 

In our recent paper [20] we have demonstrated that the simple model, originally 
proposed in Ref. jl4j, works precisely like such an escapement device. This model consists 
of a classical particle in a one- dimensional two-wave potential composed of two periodic 
in space, identical time-independent components, one of them being translated with 
respect to the other by some external force. It was found in Ref. jHj that for a constant 
or periodic in time external driving force, the particle can perform directed motion with 
a constant velocity. In our previous work j2E] we have realized that such unidirectional 
drift motion persists also in situations when the driving force is random and averages 
out to zero. We have shown that when the direction of the driving force fiuctuates 
randomly in time, upon reaching a certain level, random fiuctuations exercise a locking 
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Figure 1. Plot of the function given by Eq. with asymmetry parameters ^ = 2/5 
(sohd hne) and ^' — 1/2 (dashed line, symmetric case). 



function creating points of irreversibility which the particle can not overpass. Repeated 
in each cycle, this process ultimately results in a saltatory drift motion with random 
pausing times. 

Here we consider the case of randomly directed force in more detail; in particular, 
we generalize our previous analysis for systems in which the averaged force may have 
negative or positive non-zero values. Focussing on the overdamped limit, we exploit the 
mapping proposed in Ref.j^H]; in which the original dynamical system was viewed as 
a biased Brownian motion on a hierarchy of disconnected intervals, and derive simple 
analytical estimates for the particle's terminal velocity. These results generalize our 
previous findings and demonstrate that under certain conditions the particle may even 
perform unidirectional saltatory drift motion against the direction of the averaged 
external force. 

2. The Model 

Consider a periodic, piecewise continuous function n(x) (see, Fig.l). Without lack of 
generality, we assume that the periodicity b of this function is equal to unity, b = 1, 
so that n(a; + 1) = Il{x) for any x, and an amplitude Ho = maxn(x) = — minn(a;). 
Within each period the function Il{x) has a saw-tooth shape: 



The parameter ^ G (0,1) determines the asymmetry of the saw-tooth, with = 1/2 
corresponding to the symmetric case. We focus here on the case C, < 1/2. 

Following the model of Ref.^l], we construct now the potential V{x,'~ft) as the 
sum V{x,'jt) = n(x) + n(x — jt), where ■jt defines the external translation. Note that 




(1) 
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since both n(x) and n(x — 74) are periodic, the potential V{x,'jt) is periodic in both 
arguments, so that V{x,'jt + 1) = ^(a^^Ti) for any 7f. 

Consider now a classical particle of mass m moving in the potential V{x,'jt)- The 
deterministic equation of motion in this potential landscape reads: 

mxt + rjXt H = 0, (2) 

OXt 

where Xt denotes the particle's trajectory. Note that energy is steadily pumped into the 
system through the external translation 7f. On the other hand, the energy is steadily 
dissipated when the damping > 0, which prevents the particle from detaching from 
the system. 

Eventually, we suppose that the translation 7^ originates from an external random 
force ft- More specifically, we stipulate that the translation 7^ obeys the Langevin 
equation, 

it = ft, (3) 

ft being a random Gaussian, delta-correlated process with moments: 

Tt = fo, JU = f^ + 2D5it-t'), (4) 

where the overbar denotes averaging over thermal histories. 



3. Particle Dynamics in the Overdamped Limit 

We will restrict ourselves to the limit of an overdamped motion, ?7/[(27r) v^mllo] ^ 1 
and suppose that initially, at t = 0, a;o = and 70 = 0, so that the particle is located at 
a potential minimum. Now, in the overdamped limit, as shown in Ref.^l], the particle 
dynamics is governed entirely by the evolution of the minima of the total potential 
V{x,'yt)- Consequently, in order to obtain the trajectory Xt and estimate the statistical 
velocity V = dxl/ dt, it suffices to study the time evolution of positions of these minima. 
This process has been discussed in Ref.^l] and we address the reader to this work for 
more details. Here we outline the main conclusions of Ref.^l]: The potential V{x,'-ft) 
possesses a set of minima and the position of each minimum changes in time as the 
translation 7^ evolves. The particle, located at t = in the first minimum simply 
follows the motion of this minimum up to a certain time when it reaches the point 
X = X. The x-points are points of instability I in the [x, 74) plane, (emerging due to the 
asymmetry of functions n(x), ( < 1/2), where the corresponding local minimum of the 
potential V{x, jt) ceases to exist. As soon as the minimum at which the particle resides 
disappears, the particle performs an irreversible jump, instantaneous on the time scale of 
the translation, to one of the two still existing neighboring minima. In the overdamped 
limit, the direction of the jump is prescribed by the slope of V{x, '-ft) at this point. The 
particle moves with the second minimum until it ceases to exist and then jump to the 
next minimum and so on. 




Figure 2. Typical realization of the process Tt- The dashed lines denote the "reflection 
points" —K + 3/2 of the trajectory, while the dash-dotted lines - the points —K + ^. 



4. Evolution of the effective translation Tf 

The topology of the set X formed by the points of instability in the (x, 7t) plane has 
been amply analysed in Ref. |T3] (see, Figs.l and 2 and the discussion below) on example 
of a particle subject to the potential in Eq.(l), one component of which experiences a 
constant translation. It has been realized that the set X comprises special points —K + ^ 
and —K + 3/2, K = 1,2, . . ., at which, depending on the direction at which each of these 
points is approached by the particle, particle trajectoriy experiences some irregularities. 

To illustrate these irregularities in particle dynamics, it is expedient to introduce 
an effective translation Tt as a function obeying formally the following Langevin-type 
equation 

tt = C (Tt) + ft, (5) 

where the moments of a random Gaussian force ft have been defined in Eq.(^, while 
£(r t) is an operator, which equals zero everywhere except of the set of special points 
Tt = -K + ^ and Tt = -K + 3/2, K = 1,2,.... 

At these special points, the action of the operator C (Tt) is as follows: When Tt 
reaches the value Tt = —K+3/2, (which we will call the "reflection point"), the operator 
C (Tt) changes Tt ^ Tt — 1, i.e. shifts instantaneously its value from Tt = —K + 3/2 to 
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Tf = —K + 1/2. On the other hand, when hits the value —K + ^, the operator C iVt) 
shifts the position of the reflection point at —K + 3/2 to —K + 1/2. 

In Fig. 2 we depict a typical realization of the process Ft, defined by Eqs.dH)) and 
(0} in case when the force averages out to zero, /o = 0. In numerical simulations of 
Ft, we discretised the process defined by Eq.(jS)) taking the lattice spacing 5V equal to 
2 X 10~^ and the characteristic jump time r has been set equal to unity. Hence, the 
diffusion coefficient D = (5F)V2r = 2 x 10"^ 

Note that the operator C (Ft) encompasses all essential physics associated with the 
ratchet effect; due to the presence of the operator C (Ft) the process Ft experiences an 
effective drift in the negative direction (for C, < 1/2); in absence of C (Ft) one should, of 
course, obtain purely diffusive behavior such that Ft = and Ft = 2Dt. 

Now, the process Ft defined by Eq.Q can be viewed from a different perspective: 
namely. Ft can be regarded formally as a biased (for /q 7^ 0) Brownian motion evolving 
in time t on a hierarchical lattice composed of a semi-infinite set of independent intervals 
of length 3/2 — ^ (see Fig. 3). The Brownian motion starts at t = at the origin of the 
first interval {K = 1) and evolves freely until it either hits the right-hand-side boundary 
Ft = 1/2 (which is called, in what follows, the "reflection point"), in which case it 
gets transferred instantaneously at position Ft = —1/2 and continues its motion on the 
interval = 1, or reaches the left-hand-side boundary Ft = — 1 + ^ - the "exit point" 
and gets irreversibly transferred to the second [K = 2) interval. In the second, and etc, 
interval the process Ft evolves according the same rules. 

The trajectory of the particle Xt, a particular realization of which we depict in Fig. 4, 
follows essentially the process Ft, i.e. Xt = Ft, except for the following two situations: 
(i) when Ft appears within the intervals [—K+l, ~K + 3/2], K = 1,2, . . ., (thick lines in 
Fig. 3), the particle position stops on the left boundary of these intervals, i.e., on a given 
interval (a given K) Xt = —K + 1 when Ft is within the interval [—K + 1, —K + 3/2]. 
Therefore, the particle seizes to move for some time, this "pausing" time being equal to 
the time which Ft spends within this interval. The overall trajectory Xt is thus highly 
saltatory. 

(ii) when Ft hits the left-hand-side boundary of the K-th interval and gets 
instantaneously transferred to the next interval, the particle makes a jump from 
Xt = —K + ^ to Xt = —K. It means that although the trajectory Ft is always 
continuous, the particle trajectory Xt changes discontinuously at these points. It is 
clear, however, that such jumps of length ^ do not contribute to the average particle 
velocity Xf. Consequently, the velocity of the effective external translation V = Ft and 
average particle velocity xt coincide, i.e. V = it = Ft- 

5. Zero external bias /q. 

We turn next to the computation of the velocity V focussing first on the case of zero 
external bias, /o = 0. Evidently, even in this simplest case the map depicted in Fig. 3 
is too complex to be solved exactly. However, some simple analytical arguments can be 
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To the next interval 



Figure 3. Time evolution of Tt viewed as a diffusive process on a hierarchy of 
disconnected intervals. 

proposed to obtain very accurate estimates for the functional form of V. To do this, 
we first note that it is only when passing from the K-th interval to the next, {K + 1)- 
th interval, (which is displaced on a unit distance on the F^-axis with respect to the 
previous one), the process Tf gains an uncompensated (negative) displacement equal 
to the overlap distance of two consecutive intervals, which equals, namely, 1/2 — ^. 
Consequently, V can be estimated asV = — (1/2 — ^)/T, where T is the mean time 
which the process Fj "spends" within a given interval. To estimate the value of T, let 
us consider more precisely dynamics of the process F^ on a given interval. 

The time T matter of fact, the mean time needed for the process Fj to reach 

diffusively for the first time the left-hand-side boundary of the interval - the "exit" point 
—K-\-^, starting from the "entry" point —K + 1 + ^ (see Fig.5). It is convenient now to 
divide the i^-th interval into three sub-intervals: h = h = 1/2—^ and I2 = l/2-|-^, Fig.5. 
Evidently, since = ^3 < h, (^ > 0), for typical realizations of external random force ft, 
the following scenario holds: the process Fj evolves randomly around the "entry" point 
and first hits the right-hand-side boundary of the i^"-th interval passing thus through 
the sub-interval li; the mean time required for the first passage of this sub-interval is 
denoted as Ti. For conventional diffusion with diffusion coefficient D this time is given 
by Ti = (1/2 — C,y/2D P?]. Further on, after hitting the right-hand-side boundary, F^ 
gets reflected to the point Tt = —{2K — l)/2 and evolves randomly around this point 
until it hits the left-hand-side boundary of the interval - the "exit" point. For symmetric 
diffusion, the time T3 needed for the first passage through the sub-interval /a is equal to 
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Figure 4. Typical realization of the effective translation Fj (black line) and the 
corresponding particle trajectory xt (gray line). The dashed lines denote the "reflection 
points" X = —K + 3/2 of the trajectory Ft, while the dash-dotted hnes - the special 
points + ^. Vertical red lines show the special points at which the particle seizes 
to move. 
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Figure 5. Time evolution of the effective external translation Tt on the K-ih interval. 
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Figure 6. The absolute value of the statistical velocity V versus ^. The solid line 
gives the analytical prediction in Eq.© and the triangles depict Monte Carlo data. 
Note that we use, for notational convenience, the log-log scale. 



Ti. Consequently, T 



Ti + T3 = (1/2 — ^y/D and the particle velocity follows 



1-2^ 

In Fig. 6 we compare our theoretical prediction in Eq.© and the results of the Monte 
Carlo simulations, which comparison shows that simple theoretical arguments presented 
above capture the essential physics underlying the dynamics of the random map depicted 
in Fig. 3. 

Now, we proceed with several comments on the result in Eq.®. 
(i) First, we note that V determined by Eq.® increases indefinitely when ^ —>■ 1/2, 
although, however, V = for C, = 1/2, which is a seemingly artificial discontinuity. 
As a matter of fact, such a behavior stems from our definition of the random process 
r^; namely, by writing our Eq.© we tacitly assume that both the effective step 6T of 
the process Fj and the characteristic jump time r are infinitesimal variables, (while the 
ratio D = (5F)^/2r is supposed to be fixed and finite). On the other hand, for any 
realistic physical system ST and r may be very small, but nonetheless are both finite. 
For finite 6r and r, the process Fj can be viewed as a hopping process on a lattice with 
spacing 6T and with the characteristic time separating successive jumps equal to r. The 
hopping probabilities are symmetric; that is, the probabilities of jumping from the site 
F to either of its neighboring sites are equal to each other. For this model, we find for 
Ti and T3 the following forms (see, e.g., Ref. |27j): 

Ti = T3 = r(l+L)L, (7) 



Saltatory drift in a randomly driven two-wave potential 



10 



where L = (1/2 — ^)/(5r is the number of elementary steps ST in the interval (1/2 — ^). 
Equation ((Tj) yields, in place of Eq.®, 

Note now that for finite ST and r the velocity V in Eq. (jSJ tends to a finite value when 
^ ^1/2 (but still we have V = for ^ = 1/2). We also remark that relaxing the 
assumption of the overdamped motion and/or adding a noise term to Eq.(0), we will 
obtain a certain rounding of the behavior of V in the vicinity of ^ = 1/2 28}. As a 
matter of fact, our numerical simulations show that V becomes a smooth bell-shaped 
function with V = for ^ = 1/2 and a maximal value close to ^ = 1/2. 

(ii) Second, we notice that the velocity defined by Eqs.® or (jHl) is a monotonously 
increasing function of the parameter ^, which defines the degree of the asymmetry of the 
external potential; that is, the velocity by absolute value is maximal for low asymmetry 
when ^ — s> 1/2 and minimal, V = 2D, for the strongest asymmetry with ^ = 0. We 
would like to comment now that such a behavior relies strongly on the definition of the 
random external force ft as a Gaussian, delta-correlated noise. As ft is the property 
influenced by external processes, it might be characterized, for instance, by essential 
correlations or be a non-local in space or in time (discontinuous) stochastic process. In 
both cases, the behavior of the particle's velocity in the presence of such an external 
force would be different of that predicted by Eqs.® and (jH))- 

Suppose, for example, that ft represents the so-called delta-correlated Levy process 
with parameter fi (see, e.g., Ref. P?] for more discussion). We note here parenthetically 
that the case fi = 2 corresponds to the usual Gaussian case, which yields conventional 
diffusive motion, while the case /x = 1 describes the so-called Cauchy process. 

Now, what basically changes when we assume that ft represents the Levy process 
is the functional form of the first passage times Ti and T2. In this case the time required 
for the first passage of an interval of length 1/2 — ^ reads Ti = T3 ~ (1/2 — C,)'^ |2Z] and 
hence, the particle velocity follows 

V ~ -(1 - 20'"^ (9) 

Consequently, we may expect that the particle's velocity will be a monotonously 
increasing function of the parameter ^ only for the Levy processes with /i > 1, i.e. 
for persistent processes (which have a finite probability for return to the origin). In 
the borderline case of yU = 1 (the Cauchy process) velocity will be independent of the 
asymmetry parameter. On the other hand, for transient processes with /i < 1, (i.e. 
having zero probability of return to the origin), which are not space-filling (fractal) and 
occupy the space in clustered or locahzed patches, one may expect that the particle's 
velocity would be a decreasing function of the asymmetry parameter ^. 



6. Non-zero external bias /q. 



In this section we generalize our previous results for the generic case of non-zero averaged 
external force. We will consider both the situations when /o is negative or positive. 
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6.1. Positive bias 

Consider first the case when the bias is oriented in the positive direction, i.e. /o > 0. The 
sahent feature of this situation is that here the asymmetry of the embedding potential 
prevents the particle of moving in the positive direction, i.e. in the direction of the 
applied field /q. Hence, in this case if a drift motion of the particle takes place, - it 
would take place against the applied field! We set out to show in what follows that this 
is actually the case: the particle does travel with a constant velocity V in the negative 
direction and such an extraordinary behavior is essentially due to fluctuations in the 
process - some (exponentially small) number of its trajectories traveling against the 
external field. 

Now, we suppose that the hopping process Tt takes place on a discrete lattice with 
a spacing 6r. Each interval in Fig. 5 is thus separated into a set of L elementary steps, 
L = (3/2 — ^)/Sr. The hopping probabilities for jumps in the positive (p) and in the 
negative (g) directions are no longer equal and obey: 

p/g = exp (/5/o(57), p + q = 1, (10) 

where /3 denotes the reciprocal temperature. 

One readily notices now that in the situation under consideration the symmetry 
Ti = T3, which existed in the case /q = 0, is broken: when passing through the sub- 
interval li the process Fj follows the field, while the first passage through the sub-interval 
I3 takes place in the situation when the field /o is directed effectively against the passage 
direction. Supposing that ST and r are both finite, we have then that here P7] : 

r(l + exp(-/5/o5F)) 1/2-e 



(1 - exp ( - (3fo6r 



5T 



and 



r 



exp ( — PfoST I , , 
(^l-exp(^-/3/o5Fjj ' 



l + exp(/?/o5F)) exp(/5/o5F 



(exp (/?/o5F) - 1) L(exp [phbv) - 1 

X ( exp (/?/o(l/2 - 0) - 1) - , (12) 

respectively. 

Note that Ti grows linearly with the sub-interval length 1/2 — ^, while T3 shows 
much stronger, exponential interval-length dependence, and hence, controls the overall 
time spent within a given interval. 

Consequently, the particle's velocity in this case attains the form 

(1/2-0 sinh^ (/5/o5r/2) 



V 



cosh (/3/o5F/2) [cosh (/?/o(l - 2^ + 5F)/2) - cosh (/3/o5F/2 



(13) 
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We note again that the remarkable feature of this results is that it predicts a drift motion 
of the particle against the applied external field /q. 

In the diffusion limit, i.e. when 6r and r are set equal to zero but the ratio (5r)^/r 
is kept fixed, we find from Eq. ()13|) the following result: 

V= (1 - ^I'^'^^fS (14) 

8smh=(,(J/„(l-20/4)' 

which reduces to the result in Eq.® when /q 0. 
6.2. Negative bias 

We turn next to the case of negative bias, /o < 0. At the first glance, this case should 
not give us any surprise - one anticipates that the particle here will always be passively 
carried by the field. It appears, however, that such a behavior may take place only for 
sufficiently low asymmetries of the saw-tooth potential. For high asymmetries ^, it may 
appear that it will be more efficient for the particle to travel through the shorter interval 
(/i. Fig. 5) against the field, than to be carried passively through the longer interval I2. 
It means that exponentially small number of trajectories travelling against the field will 
still matter, as in the case of positive bias. Consequently, this will give rise to a cusp in 
the overall dependence V{^). 

We turn now to the analysis of the characteristic passage times involved in this 
case. We start first with the passage through the sub-interval ^3, in which a particle 
appearing at some time at the point —{2K — l)/2 should travel to the "exit" point 
—K + ^. The passage through this sub-interval is favored by the field /o and the first 
passage time is given by 

-(l+exp (Wr))^i/2-e 

(15) 



(1 - exp [pfo6T 



Now, the question what is the mean time typically needed to reach this point 
— {2K — l)/2 starting from the "entry" point —K -|- 1 -|- ^ is a bit more delicate: it 
can be reach either by the field-enhanced diffusion through the larger sub-interval I2 or 
by diffusion against the field /o through the smaller sub-interval /i to the right-hand-side 
boundary of the interval (the "reflection" point) followed by an instantaneous transfer 
to the point —{2K — l)/2 (see. Fig. 3). Two respective mean first passage times are 
given here by 



_r(l + exp(/3/o5rjJ^l/2 + e 



(1 - exp {pfosr 



6T 
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and 



^ _ r(l + exp {pfoSr)) (exp (/?/o(l/2 - Q)) - 1 _ 1/2-^. 
' " (l - exp {(3fo6T) ) ^ (l - exp ( - /3/o5r) ) ST \ 



The particle velocity V thus can be estimated as: 

V.--JViz« (18) 

Ts + rmn{Ti,T2) 

If Ti < T2, we would evidently recover our previous result in Eq. lfT^ up to the reversal 
of the sign, /q — /o- In the opposite case, Ti > T2, one finds 



^, {1/2 -06T r 2(5rexp(/?Mr) 

r coth ( - f3fo6T/2) ^ 1 - exp (/3/o5r) 



x(l- exp (/3/o5r/2) cosh ((/3/o5re))] ' (19) 

Now, on comparing Ti and T2 we find that these two mean first passage times are 
equal to each other when the following equality holds: 

Hence, for fixed PfoST, Ti < T2 and the particle velocity is described by Eq. lfT^ with 
/o changed for — /o when 

e > (21) 

Otherwise, T2 < Ti and the particle velocity obeys Eq. p9|) . This means that the 
overall dependence of the particle velocity on the asymmetry parameter will have a 
cusp at ^ = C,crit- Note also that ^crit is a monotonously increasing function of the 
parameter (— /3/o5r), which equals zero when PfoST = and rapidly approaches 1/2 
when {—/SfoSr) 00. Hence, ^ should be sufficiently close to its maximal value 1/2 in 
order that the behavior predicted by Eq.lfT^ takes place in the case of a negative bias. 



7. Conclusion 



In conclusion, we have studied dynamics of a classical particle in a ID potential, 
composed of two periodic, asymmetric saw-tooth components, one of which is driven 
by an external random force. Concentrating on the overdamped limit, we presented 
analytical estimates for the particle's velocity. We have shown that in such a system the 
particle will perform a unidirectional saltatory drift regardless of the fact whether the 
average force is zero, or directed against the particle motion. In case when the average 
force is negative and its direction coincides with the direction of the particle motion, we 
predicted a cusp-like dependence of the particle velocity on the asymmetry parameter 
of the saw-tooth potential. We have demonstrated that the physical mechanism 
underlying such a behavior resembles the work of the so-called escapement device, used 
by watchmakers to convert the raw power of the driving force into uniform impulses. 
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Here, indeed, upon reaching certain levels, random forces lock the particle's motion 
creating points of irreversibility, so that the particle gets uncompensated displacements. 
Repeated (randomly) in each cycle, this process ultimately results in a saltatory drift 
with random pausing times. Our analytical results for systems in which the random 
force averages out to zero are in a very good agreement with Monte Carlo data. 

8. Acknowledgments 

The authors gratefully acknowledge helpful discussions with M.Porto. GO thanks the 
AvH Foundation for the financial support via the Bessel Research Award. 

References 

[1] Maxwell J C 1872 Theory of Heat (Longmans, Green and Co, London) 
[2] von Smoluchowski M 1912 Phys. Z 13 p 1069 

[3] Feynman R P, Leighton R B and Sands M 1963 The Feynman Lectures on Physics (Addison- Wesley, 
Reading) 

[4] Curie M P 1894 J. Phys. (Paris) III 3 p 393 

[5] Ajdari A and Prost J 1993 C.R. Acad. Sci. Paris II 315 p 1635 

[6] Magnasco M O 1993 Phys. Rev. Lett. 71 p 1477 

[7] Astumian R D and Bier M 1994 Phys. Rev. Lett. 72 p 1766 

[8] Duke T, Holy T E and Leibler S 1995 Phys. Rev. Lett. 74 p 330 

[9] Astumian R D 1997 Science 276 p 917 
[10] Dialynas T E, Lindenberg K and Tsironis G P 1997 Phys. Rev. E 56 p 3976 
[11] Sokolov I M 1998 Europhys. Lett. 44 p 278; Phys. Rev. E 60 p 4946 
[12] Fisher M E and A.B.Kolomeisky A B 1999 Proc. Natl. Acad. Sci. USA 96 p 6597 
[13] Lipowski R 2000 Phys. Rev. Lett. 85 p 4401 

[14] Porto M, Urbakh M and Klafter J 2000 Phys. Rev. Lett. 85 p 491 

[15] Popescu M N, Arizmendi C M, Salas-Brito A L and Family F 2000 Phys. Rev. Lett. 85 p 3321 
[16] Lipowski R, Klumpp S and Nieuwenhuizen T M 2001 Phys. Rev. Lett. 87 p 108101 
[17] Rousselet J, Salome L, Ajdari A and Prost J 1994 Nature 370 p 446 

[18] Faucheux L P, Bourdieu L S, Kaplan P D and Libchaber A J 1995 Phys. Rev. Lett. 74 p 1504 

[19] Gorre L, loannidis E and Silberzan P 1996 Europhys. Lett. 33 p 267 

[20] Kettner C, Reimann P, Hanggi P and Miiller F 2000 Phys. Rev. E 61 p 312 

[21] Jiilicher F, Ajdari A and Prost J 1997 Rev. Mod. Phys. 69 p 126 

[22] Reimann P and Hanggi P 2002 Appl. Phys. A 75 p 169 

[23] Reimann P 2002 Phys. Repts. 361 p 57 

[24] Astumian R D and Hanggi P 2002 Physics Today (November 22) p 33 

[25] |http:// www, horologia. co.uk/ escapements. html| 

[26] Oshanin G, Klafter J and Urbakh M 2004 Europhys. Lett. 68 p 26 

[27] Lindenberg K and West B J 1986 J. Stat. Phys. 42 p 201 

[28] Oshanin G, Klafter J and Urbakh M, in preparation 



